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Parity-violating M0ller scattering measurements are a powerful probe of new physics effects, 
and the upcoming high-precision experiments will require a new level of accuracy for electroweak 
radiative corrections (EWC). First, we perform the updated calculations of one-loop EWC for M0ller 
scattering asymmetry using two different approaches: semi-automatic, precise, with FeynArts and 
FormCalc as base languages, and "by hand", with reasonable approximations. In addition, we 
provide a tuned comparison between the one-loop results obtained in two different renormalization 
schemes: on-shell and constrained differential renormalization. As the last step, we discuss the 
two-loop EWC induced by squaring one-loop diagrams, and show that the significant size of this 
partial correction indicates a need for a complete study of the two-loop EWC in order to meet the 
precision goals of future experiments. 



PACS numbers: 12.15.Lk, 13.88.+e, 25.30.Bf 



I. INTRODUCTION 



M0ller scattering is a very clean process with well-known kinematics and extremely suppressed 
backgrounds, and any inconsistency with the Standard Model will signal new physics. The next- 
generation experiment to study electron-electron scattering, MOLLER planned at JLab fol- 
lowing the 11 GeV upgrade, will offer a new level of sensitivity and measure the parity-violating 
asymmetry in the scattering of longitudinally polarized electrons off an unpolarized target to a 
precision of 0.73 ppb, and allow a determination of the weak mixing angle with an uncertainty 
of about 0.1%, a factor of five improvement over the measurement by E-158 Q. Obviously, be- 
fore we can extract reliable information from the experimental data, it is necessary to take into 
account EWC. EWC to the parity-violating (PV) M0ller scattering asymmetry were addressed in 
the literature earlier and were shown to be large A more detailed literature review can be 

found in our first work on the topic. In [Bj], we calculated a full gauge- invariant set of the 
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one-loop EWC and found the total correction to be close to —70%, with no significant theoretical 
uncertainty coming from the hadronic contributions to the vacuum polarization or other uncer- 
tain input parameters. Since it is possible that a much larger theoretical uncertainty may come 
from two-loop corrections, we investigated the importance of two- loop contribution in , by com- 
paring corrections calculated in two different renormalization schemes (RS) - on-shell (OS) and 
constrained differential renormalization (CDR) [3] - and found a difference of about 11%. This 
means that the two-loop EWC may be larger than previously thought and cannot be dismissed, 
especially in the light of precision promised by MOLLER. We divide the two-loop EWC into two 
classes: the Q-part induced by quadratic one-loop amplitudes, and the T-part which includes the 
interference of Born and two- loop diagrams. In 0, we calculated the Q-part exactly and found 
that it can reach 4%. Here, we provide a brief review of our calculations done at the one-loop 
level 0, show details of the comparison between the corrections evaluated in the OS and CDR 
schemes 0, and outline some of our calculations of higher order corrections. 



II. BORN AND ONE-LOOP CORRECTIONS 



The asymmetry between left/right longitudinally polarized electrons can be constructed in the 
following way: 

. &LL + &LR — &RL — &RR & LL ~ &RR m ^ 

Alr. 



&LL + &LR + &RL + &RR &LL + 2&LR + &RR 

enhancing the contributions induced by PV electroweak interactions. The term a = , da a stands 

° J a cos 

for the differential cross section defined in the center of mass reference frame of incoming electrons. 
At the Born level (leading order (LO)), the asymmetry is 



2/(1 - 2/) l-4a 
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Alr = n _ 2 , , ,r 2 w , y = - , (2) 
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where s?y = sin 2 6y/ = 1 — ~ 0.24. As one can see from Eq. [2J the asymmetry is highly 
sensitive to 9yy so any deviation from the SM value will signal new physics. Obviously, before we 
can extract reliable information from the experimental data, it is necessary to include EWC. The 
cross section including one-loop matrix elements is: 



r3 ^3 



a = |-|M + Mil 2 = ^-[M Q Ml + 2ReM 1 M t + MiAf}) = a + a x + a Q (3) 

where <j\ = af SE + a\ er + af ox oc a 3 is an interference term between the Born and one-loop 
amplitudes (NLO), and the cross section ctq oc a 4 is a quadratic term of the same order as the 
two-loop contribution (NNLO). To make sure that our calculations at the one-loop level are error- 
free, we evaluate EWC using two different methods. Our first method, "by hand", is to derive 
the compact analytic expressions for the leading one-loop correction manually using appropriate 
approximations for y/s < 30 GeV and y^s > 500 GeV [5|. Our second method, semi-automated, 
is to consider a full set of graphs with no approximations using computer-based algebra packages 
jl, [Io[ and 11 1 . To make sure that we calculate a gauge-invariant set of graphs, we use two sets 



of renormalization conditions (RC): the RC by Hollik (HRC) introduced in [12| for our "by hand" 
approach, and the RC proposed by Denner (DRC) in [l3| for our semi-automated method. The 
infrared divergences (IR) are treated by the soft and hard-photon bremsstrahlung (see |5|). We 
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choose our input parameters to be the fine structure constant (a = 1/137.03599), the mass of 
the W boson {m w = 80.398 GeV) and the mass of the Z boson {m z = 91.1876 GeV). A relative 
correction to the PV asymmetry is defined as 5% = {A ( l R — A° LR ) / A° LR , with the superscript in 
5j[ corresponding to the various contributions: "weak" indicates no IR-divergent graphs, and 
"QED" indicates only IR divergent graphs treated by bremsstrahlung contribution. In order to 
see how our results compare to the literature we compare the ^ eafc ([5]) = —0.2790 for the 
y^s = 100 GeV and 5^ eafc ([4]) = —0.2787 using the same input parameters as in jij and obtain 
excellent agreement. A comparison of results evaluated by two methods can be seen on Fig. 
Hl^a). To establish if the higher-order (NNLO) contributions in a given RS are important (see 
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FIG. 1: a) Weak (solid) and QED (dashed) (lj = 0.05- \J s) corrections to the PV asymmetry in M0ller scattering at 9 = 90° 
(left plot), b) Total correction to the cross section in the OS and CDR schemes (middle plot), c) Correction to the PV 
asymmetry in the OS and CDR schemes (right plot) 

[iH]). we compare results in two RS: OS and CDR. Fig. H^b) shows the total correction to the 
unpolarized cross section 5 tot = (a tot — cr°) /a calculated in the OS and CDR schemes. In the 
low-energy region, the correction to the cross section is dominated by the QED contribution, and 
the difference between the two schemes is almost constant and rather small (~ 0.01), but it grows 
at y/s > m z as the weak correction becomes comparable to QED. As a result (see Fig. []Jc)), the 
difference between the OS and CDR corrections to the PV asymmetry can reach as much as 10%, 
so contributions from two-loop corrections could become important. 




III. TWO-LOOPS CORRECTIONS: Q-PART 

The higher-order corrections (oc a 4 ) to the electroweak Born cross section can be divided into 
two classes, Q-part and T-part. The Q-part is induced by the quadratic one- loop amplitude 
(~ MiM\) (third term in Eq. |3]) and the T-part is an interference term between the Born and 
two-loop amplitudes: ot = ^-ReM 2 MQ oc a 4 (Fig. [2]). The T-part still needs to be evaluated in 
the future, but we can provide some results for the Q-part in this paper. The cross section for the 
Q-part can be divided into two terms: uq = <7q + <Jq. The first term, Oq, is an IR-divergent and 

regularized part of the cross section and the second term, Oq = (^ ) 2 5{ ■ o"o, is a finite contribution. 
The IR-divergent part has the following structure: 

^ = ^M 1 At (M 1 A + 2M/) = lg) 2 Re[^(5 1 A + 25f)] ■ a , (4) 



FIG. 2: Representative two-loop graphs for the M0ller scattering. 
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FIG. 3: Bremsstrahlung treatment of IR divergences in the Q-part. The top plot (a) represents interference between emission 
from one-loop (shaded bubble) and Born graphs. The bottom plot (b) is the two-photon emission amplitude squared. 



where 



4 In 4= 



1 +in\. Since the Q-part contains terms of order oc In 2 ^= it deserves 



a special attention. To treat the IR divergences, we have to account not only for photon emission 
from one- loop diagrams but also include a complete treatment of the two-photon emission (Fig. [3]). 
A half of the bremsstrahlung contribution in Fig. [3(a) and (b) is responsible for the treatment 
of IR divergence in the Q-part and the other half for the T-part. We take the maximum energy 
of the emitted soft photon to be u = 0.05 • y/s. The bremsstrahlung cross section for Q-part is 
derived from the Fig. [3(a) as: 
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Here, <r 7 is the total photon emission cross section and Oq is the one-photon bremsstrahlung term 
from the Q-part. The two-photons emission for the Q-part (oq 7 ) is derived from Fig. [3(b): 
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Combining Eqs. HJ [5] and [6] gives the final result for <Jq + <Jq + <Jq free from nonphysical parameters 
with the regularization parameter A cancelled analytically. Detailed calculations can be found in 
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FIG. 4: The left plot (a) shows the angular dependence of the correction induced by the quadratic part only (A a) and the 
right plot (b) shows the energy dependence of the "one-loop" , "one-loop + Q-part" and quadratic corrections. 



As one can see from Fig. H](a), the correction induced by the Q-part (A^ = {A l L ^° op+Q — 
A 1 lr OOP ) / A° LR ) can reach as much as ~ 4% at 90°. The energy dependence (Fig. H](b)) is nearly 
constant for i/s < rriz but increases rapidly after weak interactions become comparable to QED. 



IV. CONCLUSION 



With the one-loop corrections now under control, it is worth considering the electroweak radia- 
tive corrections at the two-loop level. One way to find some indication of the size of higher-order 
contributions is to compare results that are expressed in terms of quantities related to different 
renormalization schemes, and our tuned comparison between the results obtained in the on-shell 
and constrained differential renormalization schemes show a difference of about 11%. Although 
an argument can be made that the two-loop corrections are suppressed by a factor of an relative 
to the one-loop corrections, we believe that they can no longer be dismissed, especially in the light 
of the 2% uncertainty to asymmetry promised by the MOLLER experiment. At the MOLLER 
kinematic conditions, the part of the quadratic correction considered here can increase the asym- 
metry up to ~ 4%. For the high-energy region ~ 2TeV, a contribution from the quadratic 
correction can reach +30%. It is impossible to say at this time if the Q-part will be enhanced or 
cancelled by other two-loop radiative corrections, but we believe that its size demands a detailed 
and consistent consideration of the T-part, which is the current task of our group. 
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